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ABSTRACT 

We discuss high-angular resolution observations of ammonia toward four candidate high-mass 
starless cores (HMSCs). The cores were identified by the Balloon-borne Large Aperture Sub- 
millimeter Telescope (BLAST) during its 2005 survey of the Vulpecula region where 60 compact 
sources were detected simultaneously at 250, 350, and 500 /xm. Four of these cores, with no IRAS- 
PSC or MSX counterparts, were mapped with the NRAO Very Large Array (VLA) and observed 
with the Effelsberg 100m telescope in the NH3(1,1) and (2,2) spectral lines. Our observations 
indicate that the four cores are cold (Tk < 16 K) and show a filamentary and/or clumpy struc- 
ture. They also show a significant velocity substructure within ~ lkms~^. The four BLAST 
cores appear to be colder and more quiescent than other previously observed HMSC candidates, 
suggesting an earlier stage of evolution. 



Subject headings: 
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submillimeter — stars: formation — ISM: clouds — ISM: molecules — radio lines: 



1. INTRODUCTION 

The importance of massive (Af > 8M0) stars 
in shaping the Galactic structure and evolution 
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is well known (jZinnecker fc Yorkd 120071 ) . In re- 
cent years, a major observational effort has been 
made to identify the earliest stages of their evo- 
lution, which are not very well constrained, yet. 
Systematic studies by various groups have uncov- 
ered several high-mass proto-stellar objects (HM- 
POs), i.e., dense gravitationally bound cores in a 
pre-ultracompact (UC) HII region phase with typ 
ical temperatures T ^ 15 — 90 K l|Molinari et al 



Beuther et al.ll2002l ) 



19961 iMolinari et al.ll2002l ISridharan et al.ll2002l 



However, these surveys were carried out on 
ffi^S'-selected objects and thus there are very few 
observations of the so-called high-mass starless 
(or pre-protostellar) core (HMSC) stage, which is 
supposed to precede the HMPO phase, and their 
physical and kinematical properties have been 



barely studied. The survey by ISridharan et al 
(|2005[ ) was also based on previous observations of 



ffl^S'-selected candidate HMPOs at 1.2 mm, thus 
quite far from the emission peak of the spectral en- 
ergy distribution (SED) of the coldest (T < 15 K) 
co res. The survey of infrared dark clouds (IRDCs) 
bv lPillai et aD (|2006[ ) showed that IRDCs are cold 
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Table 1 
VLA Observing Parameters 



BLAST ID 


Source Name 


RA[J2000] 


DEC[J2000] 


Synthesized beam 


Sensitivity 










arcsec 


mjy beam~^ 


VIO 


BLAST J194106+235513 


19:41:06.5 


+23:55:13.5 


3.3 X 2.9 


4.7 


Vll 


BLAST J194136+232325 


19:41:36.3 


+23:23:24.9 


3.3 X 2.8 


4.3 


V27 


BLAST J194306+230125 


19:43:06.5 


+23:01:25.5 


3.4 X 2.7 


4.4 


V33 


BLAST J194319+232639 


19:43:19.0 


+23:26:39.6 


3.4 X 3.0 


5.9 



Note. — The coordinates represent the BLAST core peak positions and VLA phase tracking center. The 
synthesized beam and sensitivity refer to the NH3(1,1) hne. 



(T < 20 K) massive (M > 100 Mq) and have 
linewidths ~ 1 — 3kms~^. However, these au- 
thors probed regions with typical sizes > 1 pc, thus 
likely to represent pre-protoclusters. iMotte et al 



(|2007n in their unbiased survey of Cygnus X found 
17 cores qualifying as good candidates for host- 
ing massive IR-quiet protostars, driving outflows 
traced by SiO emission, but failed to discover the 
high-mass analogs of pre-stellar dense cores. 

Detailed s tudies of individual objects (e.g., 
G28.34+0.06. IWang et al.lbOOSilZhang et al.ll2009l : 
IRAS 05345+ 3157,'Fo ntani et al.l l2009t Cvgnus X, 



Bontemps et a l. 2009) are now contributing to set 
further constraints on the fragmentation inside 
massive dense clumps, leading to individual col- 
lapsing protostars. These works have tentatively 
identified cores in very early evolutionary stages, 
which make them good targets to study the initial 
fragmentation phases in molecular clumps. In this 
paper we describe four massive cores that appear 
to be even colder and more quiescent than similar 
objects observed in previous works, and may thus 
further constitute excellent targets for follow-up 
studies on the initial fragmentation phases. 

The Balloon-borne Large Aperture Submil- 
limeter Telescope (BLAST) has recently identi- 
fied a new and unique sample of massive, cold 
dust clumps, with characteristics sizes < 0.4 pc. 
BLAST is a 2-m stratospheric balloon tele- 
scope that observes simultaneously at 250, 350, 
and 500 ^,m. using bolometric imaging arrays 



unique in its ability to detect and characterize cold 
dust emission from both starless and protostellar 
sources, constraining the temperatures of objects 
with T < 25 K using its three-band photometry 
near the peak of the cold core SED. During the 
first BLAST science flight (BLAST05), BLAST 
conducted the first sensitive large-scale Galactic 
Plane surveys at these wavelengths. 

One of the regions observed by BLAST05 cov- 
ered 4deg^ near the open cluster NGC 6823 
in the constellation Vulpecula = 59°), at 
a distance of abo ut 2.3 kpc (see discussion in 
Chapin et al.l l2008l ) . In this region, 60 compact 



sources (< 60" diameter) were detected simulta- 
neously in all three bands. Their SEDs were con- 
strained through BLAST, IRAS, Spitzer MIPS 
and MSX photometry, with inferred dust tem- 
peratures spanning 12-40 K assuming a dust 
emissivit y index /? = 1.5, and total masses ~ 15- 



700 M0 (IChapin et al.l l2008h . At least 30% of 



these cores are new, with no IRAS-PSC or MSX 
associations. Even then, most of those with IRAS 
identifications have not been studied in detail. 
The sources detected in the BLAST bands, in 
particular those without /i?^ 5 counterparts, indi- 
cate the presence of significant quantities of cool 
dust. Thus, the BLAST observations resulted in 
an unique and uniform sample, which may contain 
a number of bona-fide HMSCs. 

Among the 60 compact BLAST sources, we se- 
lected the four coldest, massive and IR-quiet cores 



(jPascale et al.l I2008D . BLAST, until the first 



Galactic results from Herschel become available, is 
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(see Section [O]) for observations with the NRACQ 
Very Large Array (VLA) in the NH3(1,1) and (2,2) 
spectral lines. The structure of the paper is thus 
the following: we describe the observations and 
discuss the main results in Section [51 Our analy- 
sis is presented in Section [3] and finally our con- 
clusions are summarized in Section |4l 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

2.1. The Sample 

Among the starless cores candidates in Vulpec- 
ula, i.e. with no IRAS-PSC or MSX counterparts, 
we selected the four coldest (Tdust < 16 K) BLAST 
sources, VIO, Vll, V27 and V33 (see Table [1]), to 
be observed with the VLA. An analysis of Spitzer 
IRAC and MIPS images reveals that only source 
Vll might have an associated protostar, as we dis- 
cuss later (Section [33]) . 

These cores are also massive (with envelope 
masses ~ 90-200 M©) and both their absolute lu- 
minosities (ipiR ~ 50 Lq) and luminosity-to-mass 
ratios (Lpm/M < I LqMq^ , see Figure 18 of 
Chapin et al.ll2008l ) suggest an early phase of evo- 



NH. 



lution. Qualitatively, this can also be seen by po- 
sitioning the four sou rces in the L — M diagram of 
iMolinari et all ( 20081 ) (see their Figure 9), where 
they would fall in the region representing the early 
accretion phase. On the same diagram, it can be 
seen that these BLAST cores are clearly separated 
from the low-mass regime. 

None of the selected sources in Vulpecula is 
listed in th e catalog of Exten d ed Gr een Objects 
(EGOs) of ICvganowski et all (|2008h . a sample 
of massive young stellar objects outflow can- 
didates _extracted__from the GLIMPSE survey 
(pBenjamin et al. 2003). Among the Vulpecula 
sources that fall in the GLIMPSE survey area 
(which include VIO, Vll, V27 and V33) only V09 
qualifies as an EGO. This is yet another indica- 
tion of the early evolutionary phase of the BLAST 
cores in our sample. 

These four Vulpecula cores were observed with 
the VLA in the NH3(1,1) and (2,2) spectral 
lines, which are low-excitation molecular lines. 



^ The National Radio Astronomy Observatory is a facility 
of the National Science Foundation operated under coop- 
erative agreement by Associated Universities, Inc. 
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Fig. 1.— VIO. Top. NH3(1,1) VLA integrated 
spectrum of VIO. Bottom. VLA map of the inte- 
grated intensity of NH3(1,1) over main and satel- 
lite components. The first contour level is at -1-2- 
(7 and the interval between adjacent levels is 1- 
(7 (RMS = 13 mJy beam"-"^ kms~^). The values 
corresponding to the contour levels are also indi- 
cated by the horizontal lines drawn in the grey 
scale wedge to the right. The ellipse in the bot- 
tom right hand corner represents the half power 
width of the synthesized beam, while the cross in- 
dicates the nominal position of the BLAST core 
(whose uncertainty is ~ 1 arcmin) and the phase 
tracking center. The pointing position of the 100- 
m observations is shown by the asterisk and the 
large solid circle (only partially visible) represents 
the half power width of the 100-m beam. The 
dot-dashed box represents the area where the line 
emission has been integrated to generate the spec- 
trum shown in the top panel. The linear scale in 
the source is shown by the vertical bar in the left 
hand side. 
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Fig. 2.— Vll. Top. NH3(1,1) (left) and (2,2) (right) VLA integrated spectra of Vll. Bottom. VLA 
maps of the NH3(1,1) (left) and (2,2) (right) lines, integrated over the main and satellite (in the case of the 
(1,1) line) components. The boldface dashed contours in the (1,1) map indicate the specific sub-regions used 
to calculate the total mass of the source (see Section . The boldface dashed contour in the (2,2) map 
indicates the region where emission has been averaged to produce the spectrum shown in Figure [5] Other 
contours and features are as in Figure [1] (with RMS — 23mJybeam~^ kms~^ and 5 mJybeam~^ kms^^ in 
the (1,1) and (2,2) maps, respectively). 
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Fig. 3. — V27. Same as Figure [2] for V27 (with RMS = ISmJybeam ^kms ^ and 5mJybeam -'^kms ^ 
in the (1,1) and (2,2) maps, respectively). 
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Fig. 4. — V33. Same as Figure[5]for V33 (with RMS — 27mJybeani^^ kms^^ in the (1,1) map), except for 
the first contour level that here is at 3-cr. In the bottom panel only the map of the integrated (1,1) emission 
is shown (see text). 
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thus appropriate to study the cold cores in the 
BLAST catalog. In addition, the NH3 molecule 
is not expe cted to be muc h depleted in pre-stellar 
cores (e.g., lAikawa etal.l 120051 ) . although its ex- 



act abundance may depend on the co re density 
(jTafalla et al.ll200a Irlower et al.lbood) . 



2.2. VLA Observations 

The VLA ammonia observations were carried 
out in June 2008, and the array was used in its 
most compact configuration (D), with baselines 
from 35m to 1km. The NH3(1,1) and (2,2) in- 
version line s at 23.694496 and 2 3.722633 GHz, re- 
spectively (|Ho fc TownesI Il983l ). were simultane- 
ously observed in the IIF mode, with a bandwidth 
of 6.25 MHz and 128 channels, corresponding to a 
velocity coverage of about 80kms~^ and a veloc- 
ity resolution of about 0.6 km s^^. The mapped 
area in each source was ~ 2.5' x 2.5', centered 
around the nominal positions of the BLAST cores 
shown with crosses in Figures [1] to H] and Fig- 
ures [H] and [71 and with a synthesized beam FWHM 
~ 3arcsec (see Table [T|). The total time on-source 
varied from ~ 30 to ~ 40min per source and per 
line. The flux density scale was established by 
observing 3C48 and the phase calibration was en- 
sured by frequent observations of the point source 
J19259-I-21064, that had a measured flux density 
at the time of the observations of 2.1 Jy. The data 
were edited and calibrated using the Astronomical 
Image Processing System (AIPS) following stan- 
dard procedures. Imaging and deconvolution was 
performed using the IMAGR task and naturally 
weighting the visibilities. The resulting observing 
parameters are listed in Table [T] 

In order to derive the rotation temperature 
(Section 13.11) we also need to convert the VLA 
spectra in units of temperature. The conver- 
sion from flux per beam to brightness tempera- 
ture measured in the synthesized beam, Tsh, is 
obtained using the relation (at the frequency of 
the NH3(1,1) line): 



T,b[K] =2.18 



_Ft,[mJy/beam] 

^min [ 1 ^max [ ] 
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Fig. 5.— NH3(2,2) VLA spectra toward Vll (top 
panel) and V27 (bottom panel) , obtained by inte- 
grating the emission within the boldface dashed 
contours shown in the (2,2) maps of Figures [2] and 

El 



where ^min and ^max indicate respectively the mi- 
nor and major axes at half power of the elliptical 
synthesized beam. 
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Fig. 6. — Vll. VLA map of NH3(1,1) main line toward Vll, in four adjacent velocity channels (shown in 
kms~^ in the top left corner of each panel). The lowest contour level is 2-(t. The insets show the synthesized 
beam at half power and the cross indicates the nominal position of the BLAST core and the phase tracking 
center. 
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2.3. Effelsberg 100-m Telescope 

Single-point spectra of the four BLAST cores 
were also obtained with the MPIfR, lOO-m tele- 
scope and its 1.3-cm primary focus receiver. The 
pointing positions used for the lOO-m telescope 
were different from those used with the VLA and 
are shown as asterisks in Figures [1] to SI Also 
shown in these figures is the contour represent- 
ing the half power width of the 100-m beam. 
The NH3(1,1) and (2,2) lines were observed si- 
multaneously within the same band, in position- 
switched mode (with a beam-throw of 200 " along 
the RA direction), with a velocity resolution of 
0.077 km s~^ and ^ 10 minutes of on-source inte- 
gration time. 

The calibration of the 100-m data followed a 
standard procedur^ where the arbitrary noise 
tube units were converted to the antenna tem- 
perature scale, and both an atmospheric opac- 
ity and elevation corrections have been applied to 
the antenna gai The main beam efficiency was 
T^jj^i, = 0.52, the FWHM of the 100-m beam was 
39 "and the resulting noise RMS in the spectra was 

100 mK, corresponding to ^ 60mJybeam^^ (at 
the native spectral resolution of 0.077 km s^^). 

For comparison with the VLA data (Sec- 
tion 12. 5|) the 100-m spectra were also converted 
to Jy using the conversion formula: 

^.[Jy]=2.65(^4=l)V.. (2) 

where 0,„b and T,„b are the main-beam FWHM 
and brightness temperature, respectively. 

2.4. VLA Maps and Spectra 

As mentioned earlier, the VLA-D has been used 
to map the NH3(1,1) and (2,2) inversion transi- 
tions. The corresponding maps and integrated 
spectra (the integration area is represented by 
the dot-dashed boxes) are shown in Figures [T] to 
m Figures |6] and [7] show channel maps of the 
NH3(1,1) main line toward sources Vll and V27. 
The NH3(1,1) line has been detected toward all 
four BLAST cores,whereas the (2,2) transition has 



been detected in Vll, V27 and V33. This is an ex- 
cellent result, also because the positional errors of 
the BLAS T sources in Vulpecu la, which were es- 
timated bv lChapin et al.l (|2008[ ) to be - 7 - 70 ", 
could make these detections more difficult. 

We note that in V33 a reliable image of the 
NH3(2,2) emission could not be produced because 
of the low signal-to-noise ratio (SNR) and ex- 
tended emission. However, integration of the emis- 
sion in a box covering most of the NH3(1,1) emis- 
sion yielded a tentative detection of the (2,2) line, 
as shown in the top-right panel of Figure S) The 
hyperfine components have been detected only in 
the (1,1) line. One can note that, in general, the 
emission appears to arise from a clumpy and fila- 
mentary structure. In Vll and V27 this morphol- 
ogy can be seen even in the weaker (2,2) transition, 
though at the 2-(t level. 

The velocity scale shown in the spectra of Fig- 
ures [1] to [5] (and also in all subsequent spectra) is 
referred to the NH3(1,1) and (2,2) main lines. We 
note the ex cellent agreement wi th the velocities es- 
timated byjChapinieLaD (HqoS) using 13CO(1-0). 
All spectral lines are clearly very narrow and the 
VLA observations do not resolve the lines (see 
later Section 13.11 and Tables [2] and [3]). In the 
case of the (1,1) maps of sources Vll and V27 
one can note in Figures [5] and [7] the sharp change 
in the brightness spatial distribution in adjacent 
velocity channels. This indicates a significant ve- 
locity substructure within ~ 1 kms~^. In Vll and 
V27 these features are unlikely to be a consequence 
of missing flux, since with the VLA-D we recover 
most of the flux measured by the 100-ni telescope 
(see Section [23]) . 

Although the morphological difference between 
these sources could be caused by both intrinsic 
(e.g., geometrical effects and/or different evolu- 
tionary phases) and observational effects, all these 
cores present very similar features. In fact, they 
all show an internal structure, with both smaller 
cores, or fragments, and an inter-core emission, 
which also appears to be "filamentary" . The frag- 
ments inside the cores have sizes < 0.05 pc and 
we are unable to resolve structures smaller than 
~ 0.034 pc (7100 AU). 



^ |http://www.mpif r-bonn.mpg.de/div/effelsberg/| 
caiibration/calib. html 

,http: / /www.mpifr-bonn.mpg.de/staff/tpillai/efF-calib/eff-calib.htmll 
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Table 2 

Velocity, Linewidth and Total Optical Depth from Method NH3(1,1). 







VLA= 








100- 


m 




Source 


T 






T 


T 

ex 






T 




K 








K 








VIO 


4.2 ± 1.2 


33.52 ±0.07 


< 1.0 


1.8 ± 1.0 


3.23 ±0.25 


33.69 ±0.02 


0.41 ±0.05 


3.0 ±0.9 


Vll 


5.0± 1.1 


28.24 ±0.04 


< 1.0 


2.0 ±0.7 


3.97 ±0.27 


28.12 ±0.01 


0.83 ±0.03 


2.1 ±0.3 


V27 


5.5 ±2.2 


27.00 ±0.03 


< 1.0 


0.7 ±0.5 


4.63 ±0.40 


26.98 ±0.01 


0.70 ±0.03 


1.6 ±0.2 


V33 


4.3 ±0.4 


22.21 ±0.03 


< 1.0 


2.5 ±0.4 


4.97 ±0.25 


22.26 ±0.01 


0.66 ±0.02 


2.2 ±0.2 



^The VLA line parameters represent average values, estimated by integrating the NH3(1,1) and (2,2) line emission 
in the dot-dashed boxes shown in Figures [1] to H) 



2.5. Effelsberg Spectra 
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Fig. 8.— VIO. NH3(1,1) spectrum toward VIO 
taken with the Effelsberg 100-m telescope. The 
(2,2) emission was not detected. 




The Effelsberg telescope was used to observe 
the four BLAST sources toward their central posi- 
tions at higher spectral resolution. The resulting 
spectra are shown in Figures |8] to 111! where we 
note that the hyperfine components of both the 
inner and outer satellites are partly resolved, an 
indication of the narrow linewidths in this source. 

For comparison with the VLA, the spectra ob- 
tained with the 100-m telescope were resampled 
at the same spectral resolution as the VLA, and 
the results are shown in Figures [12] to [151 The 
flux measured by the VLA is almost always (ex- 
ceptions are the NH3(2,2) spectra toward Vll and 
V33) lower than that observed with the 100-m 
telescope, probably as a result of both low SNR 
in the VLA spectra and extended emission that is 
filtered out by the interferometer. If we compare 
the main-line emission, the amount of the single- 
dish flux lost by the VLA varies from ~ 20 — 30 % 
in V27 and Vll to ~ 50 % in VIO and V33. In Vll 
and V33 the VLA recovers most of the NH3(2,2) 
single-dish flux. 

3. ANALYSIS 

3.1. Derivation of Physical Parameters 
with METHOD NH3(1,1) 

The NH3(1,1) and (2,2) inversion lines show 
an electric qu adrupole and magn etic hyperfine 
structure (see IHo fc Towne^ Il983 for a review). 
The (1,1) spectra were fitted using a non-linear 



10 



E 0.5 - 




10 20 30 40 50 



V, 



[km 



Fig. 9.— Vll. Top. NH3(1,1) spectrum to- 
ward VI 1 taken with the EfFclsberg lOO-m tele- 
scope. Bottom. NH3(2,2) spectrum toward Vll 
smoothed at a velocity resolution of 0.3 km s^^. 
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Fig. 10.— V27. Same as Figure m for V27. 
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Fig. 11.— V33. Same as Figure m for V33. 



least-square method (METHOD NH3(1,1) of the 
CLASS program) which takes into account the 18 
hyperfine components. This method can deter- 
mine the optical depths and hnewidths assuming 
that all components have equal excitation temper- 
atures, that the line separation is fixed at the labo- 
ratory value and that the linewidths are identical. 

The Nll3(2,2) lines are weaker and the sensitiv- 
ity of these observations is less than the sensitiv- 
ity of the NH3(1,1) data. Thus, the quadrupole 
hyperfine structure of the Nll3(2,2) transition is 
not visible when the emission is averaged over the 
same large area as the (1,1) line, shown by the 
dot-dashed boxes in Figures [1] to IH The hyper- 
fine structure of the (2,2) line becomes visible in 
sources Vll and V27 (see Figure[5]) only when the 
emission is averaged on smaller sub-regions, such 
as the ones shown by the boldface dashed con- 
tours in Figures [2] and [3l Therefore, for the pur- 
pose of using METHOD NH3(1,1) of the CLASS 
program we always fitted the (2,2) profiles (ob- 
tained by integrating the VLA maps in the same 
region as the (1,1) transition) by simple Gaussians. 
Due to the hidden magnetic hyperfine structure, 
this tends to overestimate the intrinsic linewidths 



(jBachiller et al.lll987D . 

The results of the NH3(1,1) fitting are listed in 
Table [21 The 100-m spectra have the advantage of 
a better SNR and spectral resolution and therefore 
the derived parameters are more accurate, though 
they do not necessarily describe the same volume 
of gas. The effect of the better spectral resolu- 
tion is clearly visible in the resulting linewidths 
(see Table [2]). As mentioned earlier the VLA does 
not resolve the lines and thus the VLA linewidths 
should be interpreted only as upper limits. On 
the other hand the 100-m results show that the 
linewidths are indeed different and very narrow. 
A temperature of Tk — 14 K corresponds to a 
NH3 therma l linewidth of AV th = 0.19 kms^^ 
(A14h = ^8 In2fcrk/(17TOH), where k is the 
Boltzmann constant and toh is the atomic hy- 
drogen mass). Therefore, non-thermal motions 
definitely give an important contribution in these 
sources. The overall linewidths, however, are 
smaller (or much smaller in the case of VIO) than 
those previous ly observed toward (w armer) HM- 



POs by, e.g., iMolinari et al.l (|1996l ). who found 
me dian linewidths 1.7 — l.Skms"^, and also 
by iMotte et all (|2007l) who found linewidths > 



12 



2.0knis-i. 

Source VI 1 has the largest hnewidth, as mea- 
sured by both (1,1) and (2,2) hues. We also 
note that VI 1 is the only source where the ve- 
locity structure shown in Figure [6] is suggesting 
the presence of either two cores at separate ve- 
locities or a velocity gradient, positive S to N, of 
~ 11 kms~^ pc~^. Neither of these can be ruled 
out on the basis of the present data. Inciden- 
tally, VI 1 is also the only source with a clear posi- 
tional association with a likely protostar (see Sec- 
tion El]). 

The results of the Gaussian fits to the NH3(2,2) 
profiles are listed in Table [31 Because the Gaus- 
sian fit does not take into account the hyperfine 
structure of the line, we cannot compare directly 
the linewidths in Tables [2] and [H We also note the 
same basic behaviour among those sources with a 
(2,2) detection: source V33 has in fact the small- 
est linewidth as measured by both (1,1) and (2,2) 
lines. However, because in the VLA spectrum the 
(2,2) line is barely two channels wide (see Fig- 
ure [TSl) in Table [3] we only list an upper limit. 

3.2. Temperature 

In Vll, V27, and V33, where the (2,2) hue 
was detected, we were also able to determine 
the rotation temperature, T12 (and thus the ki- 
netic temperature T^), and the column de nsity 
using th e meth o d of Ungerechts et al. I (Il986l) and 
iBachiller et al. ( 198?! )^ In this method the re- 
quired data were: (i) the product r(Tcx — Tbg) 
for the (1,1) line, where Tex and Tbg = 2.725 K 
are the excitation and background temperatures, 
respectively, and t is the optical depth; (ii) the 
linewidth Ay(l, 1); and (in) the (2,2) integrated 
intensity. Once T12 has been estimated, the kinetic 
temperature can be determined using the analyt- 
ical expression of iTafalla et al.l (|2004l ): 



Ti: 



1 



l.lexp(-^)] 



(3) 



which is an empirical expression that fits the 
T12 — Tk relation obtained using a radiative trans- 
fer model. 

The results are shown in Table HI The VLA 
physical parameters represent average values, 
estimated by first integrating the NH3(1,1) in 
the smaller sub-regions enclosed by the boldface 



dashed contours shown in Figures [2] to |4] (see also 
Sect. I3.3p : then, the parameters obtained in each 
of these sub-regions have been averaged together 
and the results listed in Table SI Because of the 
coarse spectral resolution of the VLA spectra, 
when evaluating the physical parameters from the 
VLA data we have replaced the original linewidths 
with the 100m linewidths. Although the 100m 
linewidths represent themselves averages over a 
larger source area compared to the smaller sub- 
regions shown in Figures [2] to IH they give a better 
estimate of the Ay(l, 1) values as compared to the 
upper limits given by the VLA linewidths. In ad- 
dition, because of the weakness of the (2,2) transi- 
tion, we could not obtain separate (2,2) integrated 
spectra for the sub-regions. As a consequence, we 
have always used the NH3(2,2) spectrum obtained 
by integrating over the entire source, as shown in 
Figures [2] to m 

The small linewidths listed in Table [5] and the 
low kinetic temperatures of Table |4] again suggest 
small internal motions and therefore very quies- 
cent conditions. The four BLAST cores appear to 
be even colder and more quiescent than the HMSC 
candidates observed by ISridharan et al.l (|2005l ). 



who found in their sample of candidate HMSCs a 
median linewidth of 1 . 5kms ~^, comparable with 
that of iMolinari et al.l (|l996l ). and a median rota- 
tion temperature of 16.9 K. The average tempera- 
ture of our cores is also somewhat lower than the 
m ean kinetic tempera ture of 15 K found in IRDCs 
bv lPillai et all (l200d) . 



3.3. Column Density and Mass 

The three fitting parameters described in Sec- 
tion 13.21 may also be used to determine the col- 
umn densities, A^(l,l) and iV(2,2), of the (1,1) 
and (2,2) lines. Then, the total NH3 column den- 
sit y, A^NH-) 1 is obtain e d from T12 using the method 



of lUngerechts et al 



( 1986[ ). The total mass in 
the VLA maps has been estimated by dividing 
each source into smaller sub-regions (shown by 
the boldface dashed contours in Figures [2] to |4]) 
and determining iVNHa in each of them. The ap- 
proximate size of each region represents a trade- 
off between the need to achieve a good SNR in the 
integrated spectra of each region, and the need 
to avoid integrating over too large an area with 
no NII3 emission. For the VLA data, the values 
of temperatures and column densities averaged in 
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Table 3 

Results of Gaussian Fits to the NH3(2,2) Lines. 







VLA 






100-m 




Source 


Tsh 


Visr 


Ay 


Jinb 








K 




km s^^ 


K 




km s^^ 


VII 


0.46 ±0.09 


27.91 ±0.13 


1.4 ±0.3 


0.16 ±0.06 


28.23 ±0.19 


1.44 ±0.56 


V27 


0.24 ±0.06 


27.13 ±0.15 


1.3 ±0.3 


0.33 ±0.06 


26.96 ±0.08 


1.23 ±0.18 


V33 


0.62 ±0.14 


22.27±0.13 


< 1.0 


0.52 ±0.09 


22.36 ±0.04 


0.72 ±0.08 



Note. — The VLA fits have been performed on the integrated spectra in the area shown by 
the dot-dashed boxes in Figures [2] to ID 



Table 4 

Temperature and Column Density 



Source 


Tl2 

K 


VLA*^ 

K 


iV(NH3) 
X 10^"* cm^^ 


Tl2 

K 


100-m 

K 


A^(NH3) 
xlO"cm-2 


Vll 


10.5 ±0.5 


11.1±0.6 


21.5 ±4.7 


11.9± 1.1 


12.8 ± 1.3 


6.6 ± 1.0 


V27 


8.7±0.5 


9.0 ±0.6 


34.1 ±8.4 


14.3 ±0.6 


16.0 ± 1.0 


3.9 ±0.4 


V33 


9.2 ±0.2 


9.5 ±0.2 


37.8 ±6.7 


12.2 ±0.3 


13.1 ±0.5 


6.5 ±0.4 



'^The VLA hne parameters represent average values, estimated by integrating the NH3(1,1) 
within the sub- regions enclosed by the dashed contours in Figures[2]to|4](see text and Figure [T6|). 
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Table 5 
Size, Mass and Density 



Source 


Diameter'* 


-^BLAST 




Mtot 








arcsec 


pc 


Mq 




Mq 


Mq 


xlO^ cm" 


xlO^cm^'^ 


VIO 


7.6 


0.08 


89 


1.5 










Vll 


12.6 


0.14 


213 


10.2 


5.0 


72.0 


2.2 


31.0 


V27 


11.6 


0.13 


105 


6.7 


9.7 


138.5 


2.3 


33.4 


V33 


18.8 


0.21 


107 


9.6 


21.0 


300.7 


2.1 


30.0 



Note. — Mtot represents the total mass of the source as estimated by summing the masses 
{Med) of the sub-regions enclosed by dashed contours shown in Figures [5] to |31 (riHs) represents 
the average density among the sub-regions in which each source has been divided. 

'^Estimated from the total area of the source within the 50% contour. 



''BLAST core masses from lChapin et al.l (|2008[ ). 



■^The first and second value of both Mtot and (riHs) correspond to the values of ^[NHs] = 10~^ 
and ^[NHa] = 7 x 10~^, respectively. The errors on Aftot are derived from the errors on Nb_2 for 
each separate region, and vary between ~ 15 and 25%. 



VIO NH3(1,1) 
0.20 [ ' 

^ 0.15- j 




10 20 30 40 50 60 
Vlsr [km s"'] 



Fig. 12.— VIO. Overlay of the NH3( 1,1) spectra 
of Effelsberg (blue line) and VLA (red line) in VIO. 



these sub-regions in sources Vll, V27 and V33 are 
hsted in Table H 

As an example of the NH3(1,1) spectra and line 
fits obtained with this procedure we show in Fig- 
ure [in] the spectra obtained toward V27 by inte- 
grating the emission in the three smaller areas, en- 
closed within dashed contours, shown in Figure [21 
As mentioned already in Section [221 for each sub- 
region we have used the same NH3(2,2) spectrum 
obtained by integrating over the entire source. In 
addition, the NH3(1,1) linewidths have been kept 
fixed, and equal to the 100m values, during the 
fits to the VLA spectra obtained using METHOD 
NH3(1,1) of CLASS, such as the ones shown Fig- 
ure [161 The total mass of each source can then be 
computed by adding the masses of the individual 
sub-regions, estimated as: 

M,d^l-3S^D^Nn,mu, (4) 

where D is the sub-region diameter, mpja is the 
mass of the H2 molecule, N-^^ is the sub-region av- 
eraged column density and 1.38 is the correction 
factor for the abundance of helium and heavier 
elements in the interstellar medium. The source 
diameter is actually an equivalent diameter cal- 
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VI 1 NH3(1 ,1 



CO 




0.10 



0.05 



CD 



20 30 40 50 
Vlsr [km s"'] 



VI 1 NH3(2,2 



X 0.00 - 



-0.05 




10 20 30 40 50 
Vlsr [km s"'] 



Fig. 13.— Vll. Top. Overlay of the NH3( 1,1) 
spectra of Effelsberg (blue line) and VL A (red line) 
in Vll. Bottom. Same as above for NH3(2,2). 



>^ 0.4 



C 
QJ 




10 20 30 40 50 



Vlsr [km 

V27 NH- 



,2) 



0.20 
0.1 5 
0.10 

CO 

I 0.05 

X 

^ 0.00 

-0.05 





10 20 30 40 50 
Vlsr [km s''] 



Fig. 14.— V27. Same as Figure m for V27. 



16 



culated as -y/rig/l-lSS, where fts represents the 
sohd angle covered by the individual sub-regions 
shown in Figures [2] to HI We have evaluated 
Med for two possible values of the ammonia abun- 
dance, corresponding to the minimum and maxi- 
mum values found in Table 3 o f lPillai et al. 1 11006), 
X[NH3] = 7x 10-^ and ^[NHa] = lO"'^. The cor- 
responding values of Med are listed separately in 
Table [5] (columns 6 and 7). The particle density 
has then been calculated as: 



"-Ho 



47r(D/2)3^mH 



(5) 



where /i = 2.33 is the mean molecular weight per 
particle and mn is the mass of the hydrogen atom. 
The two values of riHj corresponding to the choice 
of XfNHa] are also listed separately in Table [5] 
(columns 8 and 9). 

From the linewidths of the observed transitions 
and the estimated source angular diameters we can 
also derive the mass required for virial equilibrium. 
Assuming the source to be spherical and homoge- 
neous (an assumption which is not, however, well 
justified in all of our sources), and neglecting con- 
tributions from magnetic field and surf ace pres; 



sure, the virial mass is given by (jMacLaren et al 
19881) : 



Mvir[M0] = 0.509 d[kpc] 6ls[arcsec] (AV^[km/s])2 

J^— 1 

where d — 2.3 kpc (jChapin et al. I l2008h is the dis- 
tance to the four BLAST sources. The estimated 
virial masses (column 5 in Table [5]) show that the 
cores are more likely to be gravitionally unstable, 
even when the largest value of ^[NHa] — 10~^ is 
used, are V33 and, to a lesser extent, V27. Source 
Vll seems to be closer to virial equilibrium, com- 
pared to V27 and V33, which would appear to be 
consistent with Vll being a proto-stellar core (see 
Section 13. 5|) . For comparison, the virial masses 
listed in Table [5] are quite smaller c ompared to the 
values estimated bv lFontani et al.l (2004) for com- 
parable sizes. However, these authors did not use 
their NH3 map and spectra to determine the virial 
mass, and they measured much larger linewidths 
in other molecular tracers. Therefore, the virial 



masses of Fontani et al.l (|2004l ) are likely to reflect 
the much larger degree of turbulence in HMPOs. 

Furthermore, apart from the various approxi- 
mations involved in the calculation of Afvir, its 



comparison with the masses obtained by the fit 
to the spectral energy distribution (with assumed 
values of the dust mass absorption coefficient at 
250 /j,m, k25Q — 10 cm^g~ ^, and a dust emissiv- 



ity index, /3 = 1.5, Chapin et al. 2008f) is further 



complicated by the fact that the BLAST05 ob- 
servations were sensitive to cold dust on a much 
larger scale (> larcmin). 

The comparison between Mtot and M^ir in Ta- 
ble [3] is complicated by the many uncertainty fac- 
tors. In addition, for a power-law density distri- 
bution of the type nB_2i'r) oc r"™, the virial mass 
obtained from Eg. (161) must be niultipl ied by a fac- 
tor f (s ee'MacLare n et al.llTgsil ). For exam- 
ple, ^bntanLgt al. (2004) found that m = 2.3 in 
the HMPO IRAS 23385-^6053 and thus their virial 
masses had to be multiplied by a factor ~ 0.35. 

3.4. Comparison with Lov^r-Mass Pre-Stellar 
Cores 

To analyze the origin of the VIO, Vll, V27 and 
V33 cores, we can compare the properties of these 
cores with the well-studied low-mass pre-stellar 
core L1544 in the Taurus molecular cloud (at a 
distance of 140 pc) that appears to be gravitation- 



ally unstable (for a review see iDotv et all l2005l 
and references therein). The LI 544 core has a 
continuum s ource with a flux d ensity of 17.4 Jy 
at 450 fj,m (IShirlev et al.l l2000l) . corresponding 
to a gas-l-dust ma ss of 3.2 M(^, and a densit y 
of 1.5 X lO^cm-^ (|Ward-Thompson et al.|[l999l) . 
The N2H+ linewidth toward LI 544 is 0.3 km s-\ 
and observations toward several low-mass starless 
cores in NH3 and N2H"'" show similar linewidths 



for both molecules ( Tafalla et al.ll2004[ ) 



In order to compare the physical properties of 
L1544 with those of the BLAST cores listed in 
Table O we n ote that the peak N H3 abundance 
determined bv lTafalla et all (|2002l ) toward L1544 
was 4.0 X 10~^. Therefore, if we consider the phys- 
ical parameters in Table [5] corresponding to the 
lower value of the NH3 abundance (7.0 x 10~^, 
close enough to that estimated in L1544), we see 
that the number density in the BLAST cores is 
quite similar to that of L1544, while their mass is a 
factor of ^ 20—100 higher, and is quite higher even 
when one considers the masses of the smaller sub- 
regions within each source, as described in Sec- 
tion [ 
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15.— V33. Same as Figure m for V33. 



T,2 = 8.2+/-0.4 K 
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Fig. 16. — VLA spectra obtained towards V27 
by selecting the three smaller sub-regions enclosed 
by the boldface dashed contours shown in Fig- 
ure El The solid line represents the fit obtained 
with METHOD NH3(1,1) of the CLASS program, 
by setting the VLA linewidths to the 100m values 
(see text). Also shown are the values of the rota- 
tion temperature (Section 13. 2|) and NH3 column 
density as obtained in each separate sub-region 
(see Section 13.31 and Table lU where the param- 
eters listed are weighted averages). 
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As far as the velocity dispersion is concerned, 
the BLAST cores have NH3 linewidths more than 
double that of L1544, with the exception of VIO 
which has a linewidth (0.41 km s~^, see Table [5]) 
only slightly higher than that of LI 544. This 
would be consistent with the BLAST cores being 
more turbulent and could also possibly hide other 
systematic internal motions (see the case of VI 1, 
Section EH). 

What would a core like L1544 appear in the 
BLAST05 map of Vulpecula? Scaling the 450 ^m 
L1544 emission to 500 /im at a distance of 2.3 kpc, 
we would expect an integrated flux density of 
about 40mJy and thus it would go und etected 



in th e BLAST05 map of Vulpecula tChapin et al 
20081 ). This is consistent with the four BLAST 



sources being separated, if they w ere plotted on 
the L — M diagram of .Molinari et al . (2008), from 
the low-mass regime, as mentioned in Section [2. II 
Table [S] shows that only in the case of VI 1 the 
total mass inferred from the BLAST observations 
is significantly larger than the mass estimate based 
on NH3 emission, given the NH3 abundance range 
that we have selected. In V27 and V33, the mass 
range obtained from our observations is consistent 
with the BLAST mass estimate. It is possible that 
in the case of VI 1 the ammonia emission results 
from only a small fraction of molecular material 
that is concentrated in a denser and more compact 
structure. A similar scenario , for e xample, has 
been found by IShepherd et al.l (|2004l ) toward the 
core G192 S3, and it could also be the case in 
V27 and V33 if the NII3 abundance were near the 
upper limit of the range listed in Table \5\ In this 
case, an interesting question would be whether the 
mass determined from the VLA NH3 observations 
represents the total reservoir from which any star 
formation will eventually take place in these cores 
or, rather, further mass will be accreted from the 
larger mass reservoir detected by BLAST. In the 
first scenario, it is unlikely that these cores will 
produce massive stars, unless the total gas and 
dust mass is much higher, as it would result if 
X[NH3] « 10-^. 

Likewise, one might ask whether the lower-mass 
fragments in these regions could at a later stage 
merge to form a more massive pre-stellar core 
that could then evolve toward a high-mass star. 
Clearly, we cannot currently determine the rela- 
tive motions of the individual fragments in the 



V 1 



23' 15 



23" 22' 45' 




37^ 35^ 
Right Ascension 



Center: R.A. 19 41 36.32 



Dec +23 23 23.5 



Fig. 17. — Three-color Spitzer image 3.6 /xm 
(blue), 8^m (green), and 24 ^m (red) of VI 1 
(at the angular resolution of the MIPS 24 /im 
map). The overlaid white, solid contours rep- 
resent the BLAST 350 /ttm emission (from 4.3 x 
10^ to 1.2 X 10^ by 1.6 X 10^ Jy srad^^). The 
black, solid contours represent the integrated 
NH3(1,1) VLA emission (from 0.039 to 0.112 by 
0.036 Jybeam"^ kms~^). 



19 



observed cores. However, we can use the NH3 
linewidths as a rough estimate of the average tur- 
bulent motion, though they do not necessarily 
represent the relative velocities of coherent sub- 
structures within the cores. If we take ~ 0.1 pc 
and O.S-s/Skms"^ as the typical fragment (de- 
projected) separation and 3D velocity dispersion, 
respectively, we obtain a time-scale ~ 8 x 10^ yr. 
Although this time-scale is comp arable to the ac- 
cretio n time-scales estimated by iBonnell fc Bate 
( 20061 ) for the formation of high- mass stars in the 
competitive accretion scenario, the most signifi- 
cant individual fragments inside the cores appear 
to be more massive than the thermal Jeans mass 
(~ IMp,), which i s incon sistent with the model 
of iBonnell fc Batd (l2006l) for competitive accre- 
tion. Alternatively, it is quite possible that we 
are observing these cores at a phase when merging 
of low-mass fragments has already occurred, and 
that the less significant fragments in Figures [T] to 
m with masses > 1 Mq , are actually the remnants 
of this accreting phase. 

3.5. Evolutionary Phase of the Cores 



As we have seen in Section [3?2] and Section [3?3l 
our cores have similar temperatures but are more 
compact, and possibly less massive (depending on 
the NH3 abundance), than IRDCs, which have a 
typical size > 1 pc. This suggests that they may 
be similar to the unresolved cores found toward 
IRDCs, but in an earlier stage of evolution com- 



pared to t he sources observed by iMolinari et al 
(|l99fil ) and lSridharanet all (120051 ). 



This conclusion is supported by the lack of an 
IRAS and MSX counterpart to these sources, as 
mentioned in Section [2Tl However, we have also 
inspected the IRAC and MIPS images and found 
that only in the case of VI 1 an IRAC/MIPS point 
source is found to be positionally coincident with 
the NH3(1,1) integrated emission. As shown in 
Figure[T7]this source has the colors expected for an 
embedded protostar. Aperture photometry con- 
firms that the SED of this infrared source is ris- 
ing between 3.6 and 24 ^m. However, we find 
no point source in the 70 /xm MIPS image. No 
other infrared source with similar colors is found 
within ~ 1 arcmin from the peak of either the 
NH3(1,1) or the BLAST emission. Together with 
the larger linewidth and a more systematic veloc- 
ity shift across the source structure, as discussed 



in Section [3TT1 these properties suggest that VI 1 
is probably the most evolved of all four observed 
cores, and is likely to harbour a protostar. 

In addition, we have also searched for possible 
radio continuum emission toward the four Vulpec- 
ula sources. In order to do this we have analyzed 
both the continuum images constructed from the 
line-free channels of our VLA observations, and 
the available maps f rom the CORNISH survey of 
the Galactic Plane (jPurceh et al.ll2008[) . We do 



not detect any continuum emission at the level of 
~ 0.4mJybeam^^, in our line-free channels maps, 
which is about the same as the noise level achieved 
at 5 GHz in the CORNISH maps. Thus, both 
the upper limit to the continuum emission (see 
Purcell et al.ll2008h and t he luminosity estini ated 
from SED fits (L ~ 50 Lq^. lChapin et al.lliooi ) im- 
ply an ionizing source much less than a B3 zero- 
age main-sequence star. The low-luminosity in the 
radio continuum could be because the central ob- 
ject has not yet developed an UC HII region, or 
because of opacity if the surroundings of the pro- 
tostar (if any) are very dense. 

4. CONCLUSIONS 

We have observed four candidate high-mass 
starless cores selected from the BLAST05 survey 
of the Vulpecula region. Our VLA-D observations 
in the NH3(1,1) and (2,2) lines suggest that these 
cores are in very early stages of evolution, prior 
to the formation of a (proto)star or proto-cluster. 
In only one of these cores, Vll, the VLA peak 
emission of NH3(1,1) is associated with an infrared 
source having the colors expected for an embedded 
protostar. The four cores are cold (Tk < 16 K), 
relatively quiescent {AV < 0.8 km s"^) but with a 
higher velocity dispersion compared to e.g., L1544, 
and show a filamentary and clumpy structure. Our 
VLA-D data have limited spectral resolution, but 
we can clearly observe a significant velocity sub- 
structure within ^ Ikms^^. 

Based on the comparison with a typical low- 
mass pre-stellar core, we find that these BLAST 
cores are more massive and more luminous. Vll 
is likely to have already formed an intermediate- 
mass (proto)star and given the mass range ob- 
tained for V27 and V33 they too could form 
intermediate- or even high-mass stars. However, a 
confirmation of this conclusion will require a more 
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accurate determination of the NH3 abundance in 
these cores. 
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